Relic gravitational waves: latest revisions and preparations for new data 
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The forthcoming release of data from the Planck mission, and possibly from the next round of 
Wilkinson Microwave Anisotropy Probe (WMAP) observations, make it necessary to revise the 
evaluations of relic gravitational waves in the existing data and, at the same time, to refine the 
assumptions and data analysis techniques in preparation for the arrival of new data. We reconsider 
with the help of the commonly used CosmoMC numerical package the previously found indications 
of relic gravitational waves in the 7-year (WMAP7) data. The CosmoMC approach reduces the 
confidence of these indications from approximately 2a level to approximately la level, but the 
indications do not disappear altogether. We critically analyze the assumptions that are currently 
used in the cosmic microwave background (CMB) data analyses and outline the strategy that should 
help avoid the oversight of relic gravitational waves in the future CMB data. In particular, it 
is important to keep away from the unwarranted assumptions about density perturbations. The 
prospects of confident detection of relic gravitational waves by the Planck satellite have worsened, 
but they are still good. It appears that more effort will be required in order to mitigate the 
foreground contamination. 



I. INTRODUCTION 

Cosmology relies on extrapolations. It is natural to 
extrapolate the picture of approximate homogeneity and 
isotropy seen in our patch of the Universe to other places 
and times. Having assumed the usual framework of a 
homogeneous and isotropic hot big bang, we do not en- 
counter any real contradictions or paradoxes requiring 
drastic solutions. But we do encounter an initial cosmo- 
logical singularity [l|. The singularity is not an ultimate 
answer, it is only a sign of limited applicability of the 
currently available theories. The limit of applicability is 
probably set by the Planck parameters. It seems logical 
to suggest that our Universe came into being as a con- 
figuration with a Planckian size and a Planckian energy 
density, and with a total energy, including gravity, equal 
to zero (see Q and references therein). 

It is here where a real problem arises. The newly cre- 
ated classical configuration cannot reach the averaged en- 
ergy density and size of the presently observed Universe, 
unless the configuration experienced a kind of primordial 
kick. During the kick the size of the newly born Uni- 
verse (for simplicity, one can think of the curvature ra- 
dius of a closed universe) should increase by many orders 
of magnitude without significant changes in the energy 
density of the whatever substance that was there. Such 
an expansion cannot be driven by normal types of matter 
which we presently know. The kick should be driven by 
something more exotic: a lucky version of a scalar field 
(inflaton) @ or something like a conformal anomaly [|| 
or something from the "theory of everything" . 

The strong variable gravitational field of the very early 
Universe inevitably generates gravitational waves and, 
under certain extra conditions, density perturbations and 
rotational perturbations [5j- The generating mechanism 
is the superadiabatic (parametric) amplification of the 



zero-point quantum oscillations of the respective degrees 
of freedom, or, in more technical terms, the Schrodinger 
evolution of the initial vacuum (ground) state of the cor- 
responding time-dependent Hamiltonian into a strongly 
squeezed vacuum (multiparticle) quantum state [f| . The 
generated gravitational waves and density perturbations 
leave observable imprints in the temperature and polar- 
ization anisotropics of the cosmic microwave background 
radiation. By studying the CMB correlation functions 
TT, TE, EE, BB we can learn about the birth of the 
Universe and the initial stage of cosmological expansion. 

The simplest assumption about the initial kick is that 
its entire duration can be described by the scale factor 
with one fixed power-law dependence [5] : 



a( V ) = l \ V \ 1+ ^ 



(1) 



where l D and /3 are constants, j3 < — 1. Then, the gener- 
ated primordial spectra of metric perturbations describ- 
ing gravitational waves (t) and density perturbations (s) 
(for more explanations, see have the universal power- 
law forms 



(2) 



P s (k) =A S [- 



where we will be using fco = 0.002Mpc _1 . According to 
the theory of quantum-mechanical generation of cosmo- 
logical perturbations [5] the spectral indices are approxi- 
mately equal, n s — 1 = n t = 2(/3 + 2), and the amplitudes 
(At) 1 / 2 and (As) 1 / 2 are of the order of magnitude of the 
ratio Hi/ Hp\, where Hi ~ c/l a is the characteristic value 
of the Hubble parameter during the kick. 

The presently available CMB and large-scale struc- 
ture observations cover a huge interval of scales, rang- 
ing from wave numbers k « 0.0002Mpc _1 and up to 
k « 0.2Mpc _1 . There is no any particularly fundamen- 
tal reason why a single power-law evolution (UJ should 
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be valid during a long interval of time and, consequently, 
exact power-law spectra @ valid in the interval of wave- 
lengths where the ends of the spectrum differ from each 
other by a factor 10 3 . We do not know such spectra in 
physics and astrophysics. In cosmology, one can also ex- 
pect some deviations from strict power-laws (J21) 0j B- 

If the continuous power spectrum can be approximated 
by two power-law pieces, one can write 



Ps(k) = A™ 



P s (k) = 4 2 > 



(2) ( 

where A\ = A 



k 

1) k±\ ° S 



k < hi, (3) 
k > h, (4) 



and 



Pt(k) 



A 



(i) 



Pt(k) = A 



(2) 



where A 
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k 
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h), 

(!) „( 2 ) 



k < fci, (5) 
k > fci, (6) 
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Obviously, we return to Eq.Q if the spectral indices 
in the two pieces are equal to each other. 

Alternatively, one can postulate a simple law of "run- 
ning" of the spectral index, 



n s (k) = n s (k ) + a s ln(fc/fc ) 



(7) 



where a s is a constant. The spectral index n s {k) — 1 
for density perturbations is defined as n s {k) — 1 = 
[d \nP s (k)/d Ink], and similar definition is used for 
nt(fc). Then the power spectrum for density perturba- 
tions takes the form 



P.(k) = A e (ko) [ - 



^ \ n s (fco) — l+jcts ln(fe/fc ) 



(8) 



and a similar formula holds for gravitational waves. In 
our calculations below we will be testing both options, 
(j3])-(jl|) and (0, of spectral deviations. 

Following the (quite unfortunate) tradition, we are 
characterizing the amount of relic gravitational waves 
(GW) by ratios involving density perturbations (DP), 
rather than by the gravitational wave amplitudes di- 
rectly. One measure is the quadrupole ratio 



R 



g£ T 2 (gw) 

cji 2 (d P y 



another - the so-called tcnsor-to-scalar ratio 

AM 
A s (k )' 



(9) 



(10) 



By definition, r and R cannot be negative. Some ele- 
ments of numerical software are better adjusted to one 



quantity than to another. We will be using both of them 
remembering that r « 2R for a quite wide class of mod- 
els 0. 

Contrary to the theory of quantum-mechanical gener- 
ation of cosmological perturbations, the inflationary the- 
ory claims that the power spectrum of density perturba- 
tions should be many orders of magnitude larger than the 
power spectrum of gravitational waves. The inflationary 
theory attempts to apply to density perturbations the 
same mechanism of superadiabatic (parametric) ampli- 
fication that was originally worked out for gravitational 
waves But the so-called classic result of inflation- 
ary theory states that the A s should be arbitrarily large 
in the limit of dc Sitter inflation (i.e. when the pertur- 
bations are generated with spectral indices n s = 1 and 
n t = 0): A s « A t /e, where e = —H/H 2 is zero for de Sit- 
ter expansion with any Hubble parameter H (i.e. for any 
"energy scale of inflation"). The existence and properties 
of primordial density perturbations are consequences of 
the quantum-mechanical generating mechanism, and not 
of the inflationary theory. The contribution of inflation- 
ary theory to this subject is the (incorrect) prediction of 
divergent A s , that is, the (incorrect) prediction of am- 
plitudes that can be arbitrarily larger than Hi/Hp\. For 
inflationary considerations about density perturbations, 
see for example and recent books on cosmology. 

Inflationists rarely suffer from the problem of "too lit- 
tle", they usually have the fine-tuning trouble because 
they generate "too much" . They have started from the 
threat to "overdose the Universe" and they continue this 
trend till now. To make the wrong theory look "consis- 
tent" , inflationists convert their prediction of arbitrarily 
large A s into the prediction of arbitrarily small A t . As a 
result, the most advanced string-based inflationary theo- 
ries predict the negligibly small amounts of gravitational 
waves, r sa 10~ 24 or less. In other words, inflationists 
attempt to pay by gravitational waves for their incorrect 
treatment of density perturbations. (For a more detailed 
criticism of inflationary theory, see the last paper in [f|.) 
Obviously, we do not accept the inflationary theory. We 
allow the use of the inflationary "consistency relation" 
r = — 8n t only in Sec [Til where we repeat the derivations 
of the WMAP Collaboration [TJ], [Tj who used this re- 
lation. 

The previous search 0] for relic gravitational waves 
in the WMAP data was based on the simplified likeli- 
hood function which, among other things, neglects the 
data and noise correlations at the lowest multipoles. 
The approaching release of data from the Planck mis- 
sion [12| and, possibly, from the next round of WMAP 
observations, makes it necessary to test by other nu- 
merical techniques the previously found indications of 
relic gravitational waves. The WMAP team and other 
groups routinely employ the CosmoMC software [TJj]. 
The data analysis constructions in Q are based on the 
genuine Wishart probability density function, whereas 
the WMAP likelihood function, as a part of CosmoMC 
package, is based on semi-Gaussian approximations to 
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this distribution. The new analysis with CosmoMC is 
certainly different from the previous one, but we are re- 
luctant to say that it is necessarily more accurate from 
the physical point of view. However, we realize that the 
CosmoMC code is the approach that most of the groups 
are pursuing now and will be using in the future with 
new data. Therefore, it is important to revise the previ- 
ous evaluations of relic gravitational waves with the help 
of CosmoMC, and to make revised forecasts for Planck 
observations. This is the main purpose of this paper. 

In Sec|n]we rederive the results of WMAP Collabora- 
tion [loj. using their assumptions and methods. In 
Sec lIIII we apply the CosmoMC package to repeat the 
previously performed search for relic GW at lower-^ mul- 
tipoles Q . We show that the indications of G W diminish 
but do not disappear altogether. In SeclIVI we demon- 
strate that the conclusions about n s and r depend on 
the interval of multipoles utilized in the data analysis. 
It is unwise to postulate constant spectral indices in a 
huge interval of accessible wavelengths and multipoles. 
In Sec|V] we explore a hypothesis of primordial spectra 
consisting of two power-law pieces. It is shown that the 
hypothesis is consistent with the data and with the re- 
sults discussed in previous Sections. SeclVIlis devoted to 
refined forecasts for the Planck mission. The prospects of 
confident detection of relic gravitational waves by Planck 
are still reasonable. 



II. REPEATING THE WMAP7 ANALYSIS 

First of all, with the help of CosmoMC, we repeat the 
WMAP7 calculations. We want to be sure that we ob- 
tain the same results when we make the same assump- 
tions. Obviously, we include gravitational waves from the 
very beginning, we do not consider the "minimal" model 
of WMAP [l0| where gravitational waves are voluntar- 
ily excluded from the outset. In Sec. Ill A[ the spectral 
indices are assumed strictly constant in the entire inter- 
val of relevant wavelengths, that is, together with the 
WMAP team, we adopt Eq.©. In Sec. UTBl we allow 
for the "running" of n s , that is, together with WMAP, 
we use Eq.([7]). Moreover, in these two subsections, and 
only there, instead of the correct relation n t = n s — 1, 
we use the incorrect "consistency relation" n t — — r/8 
for reducing the number of unknown parameters. This is 
done exclusively because the WMAP team has done this, 
and we want to rederive their results. Since r cannot be 
negative, the inflationary theory does not allow positive 
n t . 



A. Constant spectral indices 

We repeat the WMAP team analysis of the 7-year 
CMB data with the help of the CosmoMC sampler. The 
space of free parameters subject to evaluation consists of 
four background parameters (flbh 2 , il c h 2 , t, 9) and three 



perturbation parameters (ln(10 10 A s ), n s , r). These are 
the standard CosmoMC parameters for a flat ACDM cos- 
mology. As mentioned, we also adopt Eq.© and atop of 
that, n t = —r/8. We call this set of assumptions the 
Case I. 

The analysis takes into account the observed CMB 
data for all four angular power spectra: TT, TE, EE 
and BB [jjj . The range of the used multipoles is 
2 < I < ^max, where Cax = 1200 for TT spectrum, 
800 for TE spectrum, and 23 for EE and BB spectra. 
The maximum likelihood values and marginalized distri- 
butions are found for all seven free parameters, but we 
are showing the results only for n s and r because we are 
mostly interested in them. 

The results of the CosmoMC calculations for the Case 
I are summarized in the upper row of Table U and in 
Fig. [T] The seven-dimensional (7d) ML values of n s and 
r are shown in the column "Maximum likelihood" . The 
one-dimensional (Id) marginalized distributions for n s 
and r are shown in Fig. Q] by olive curves, marked also by 
symbol (1). (The maximum of Id likelihoods is always 
normalized to 1.) The obtained results are very close to 
the WMAP findings [10|. We see that the assumptions 
of the Case I lead to the conclusion that the primordial 
spectrum is red, n s < 1. Although the uncertainties are 
still large, the n s < 1 is the preferred outcome of the 
WMAP analysis and of our repetition here under 
the same assumptions. The n s < 1 would mean that 
the function P s {k) is infra-red divergent at very small 
wavenumbers k. The maximum likelihood (ML) values 
of r, both, for the 7d and Id posterior distributions, are 
very close to zero. The 95% confidence limit r « 0.38 
is usually interpreted as the upper limit for the possible 
amount of relic gravitational waves. 



B. Running spectral index 

Here, together with the WMAP Collaboration, we in- 
crease the number of free parameters from seven to eight 
by adding the parameter a s and adopting Eq. (J?} . We call 
this set of assumptions the Case II. Other parameters, 
observational data, and CosmoMC numerical techniques 
are exactly the same as in the Case I. The parameter a s 
does not have to be a nonzero number. It is the maxi- 
mum likelihood analysis of the data that determines its 
value. 

The results of the analysis for the Case II are shown 
in the lower row of Table U and in Fig. [T] The eight- 
dimensional (8d) ML values of n s (ko), a s and r are shown 
in the column "Maximum likelihood" . The posterior Id 
distributions of n s (ko) and r are plotted by magenta 
curves, marked also by symbol (2). 

Together with the WMAP team, we see that the pri- 
mordial spectrum in the Case II becomes blue at long 
wavelengths, that is, the 8d and Id ML values of n s (ko) 
are larger than 1. The power spectrum P s (k) is no longer 
divergent at very small k. Since the ML a s is negative, 
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the spectrum gradually turns over from the blue at long 
wavelengths to red, n s (k) < 1, at shorter wavelengths, 
see also Eq.Q. There is a considerable decrease of n s (k) 
on the way from the long- wavelength end of the spectrum 
at k rs 0.0002Mpc _1 to the short-wavelength end of the 
spectrum at k « 0.2Mpc _1 . (We restrain from a discus- 
sion of disastrous consequences that the values n s (k) > 1 
bring forth to the inflationary theory.) 

In comparison with the Case I, the Id distribution of 
r becomes significantly broader, and it is almost flat at 
small values of r, where r G (0, 0.12). The 8d ML value 
of r is r — 0.06. This value is 7.5 times larger than that in 
the Case I, but one should take this number with caution. 
When the distribution is almost flat, the ML point can 
be accidental. It is clear from the shape of the magenta 
curve in Fig. [T] that the assumptions of the Case II are 
consistent with the hypothesis of no gravitational waves, 
r = 0. Nevertheless, the probability of large values of 
r has increased, as compared with the Case I. The 95% 
C.L. has risen up to r » 0.43. 



III. EVALUATIONS OF RELIC 
GRAVITATIONAL WAVES FROM LOWERS 
CMB DATA 

The purpose of this Section is to revise by the Cos- 
moMC code the previously found indications of relic 




in the WMAP 7 data. It was 
15j that the relic GW compete 



gravitational waves 
stressed many times 
with density perturbations only at relatively small mul- 
tipoles £. It is dangerous to include high-£ CMB data in 
the search for gravitational waves, as the spectral indices 
may not be constant. Conclusions about the amount of 
relic gw in the interval 2 < £ < 100 depend on the contri- 
bution of density perturbations to these multipoles, but 
the assessment about the participating A s and n s may 
be wrong, if it is built on the hypothesis of constant spec- 
tral indices in a huge interval of the observed multipoles. 
The relatively small number of data points in the interval 
2 < £ < 100 brings forth the increased uncertainty in the 
evaluation of relic GW, but it seems to be wiser to live 
with larger uncertainty (and wait for better data) than 
with artificial certainty based on wrong assumptions. 

The arguments put forward in @ , as well as the results 
of Sec lIIBl hint at possible deviations of the underlying 
spectra of primordial perturbations from strict power- 
laws. We perform a special investigation of this issue in 
SecED 

No doubt, in a search for relic GW, one should always 
check for the presence of residual systematic effects and 
alternative explanations. The lower-f multipoles attract 
attention for a variety of reasons, including purely in- 
strumental deficiencies, and they all should be examined 
(see, for example, [H, 0, 0). The helpful signatures 
of gravitational waves are a nonzero B-mode of CMB 
polarization [l!| and a negative TE cross-correlation at 
lower multipoles [20l ]. 



In Sec. Ill A we make exactly the same assumptions 
as in @, including the fixed best- fit values of the back- 
ground parameters. But the new results for the perturba- 
tion parameters A s , n s and r follow from the CosmoMC 
technique rather than from the procedure of Ref. [6] . We 
generalize this search in Sec lIII Bl where the background 
parameters are not fixed but are derived together with 
the perturbation parameters. For this purpose we use 
the CosmoMC facility allowing to combine external data, 
not affecting the perturbation parameters, with the CMB 
data at 2 < £ < 100. 



A. Fixed background parameters 

To work exactly with the same assumptions as in 
Ref.[6j, we fix the background parameters at their best- 
fit values of ACDM cosmology fl b h 2 = 0.02260, 
fl c h 2 = 0.1123, Q A = 0.728, r = 0.087. The Hubble 
constant h = 0.704 is a derived parameter. The free pa- 
rameters subject to evaluation by the CosmoMC code 
are ln(10 10 A s ), n s and r. The spectral indices are re- 
lated by rit = n s — 1. The CMB data are used only up to 
£max — 100. We call this set of assumptions the Case III. 

There exists some awkwardness in our choice of the 
background parameters, as these are t he p arameters that 
were derived by the WMAP team [l0| under the as- 
sumptions of complete absence of gravitational waves 
and strictly constant n s in the entire interval of partic- 
ipating wavelengths. We know from our previous expe- 
rience [6j that the background parameters, if changed 
not too much, do not greatly affect the results for the 
perturbation parameters. Nevertheless, for safety, we 
explicitly explore the issue of background parameters, 
and allow them to vary, in the more general approach of 
Sec lIII Bl This will remove the aforementioned awkward- 
ness in choosing and fixing the background parameters. 

By running the CosmoMC code, we find the following 
3d ML values of the perturbation parameters r, n s and 
lntlO 10 ^): 



0.285, 



1.052, ln(10 lu A s ) = 3.023, (11) 



and Tit — 0.052. We also derived the marginalized Id 
results for these parameters (see also Table HI)) . 



n s = 1.064 



+0.058 
-0.059' 



ln(10 iu A s ) = 2.996t°;^. (12) 



Following the convention of (Tpj , we quote here and below 
the mean values of the Id likelihood functions, and the 
uncertainties refer to the 68% confidence intervals. The 
Id likelihood functions for n s and r are plotted by red 
lines, also marked by (3), in the upper and lower panels 
of Fig. [TJ respectively. 

The red curve in the lower panel shows a clearly visi- 
ble broad maximum at r = 0.2. Unfortunately, the peak 
point r = 0.2 is not strongly separated from r = 0, it is 
nearly (slightly less than) la away from r = 0. Speaking 
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TABLE I: Results for n s , a a and r in Case I and Case II 





Maximum likelihood 


1-d likelihood 




n s 




r 


n s 




r (95%C.L.) 


Case I 


0.967 




0.008 


u.y»±_ 020 




r < 0.379 


Case II 


1.061 


-0.043 


0.060 


l.UDO_ 05g 


-0-039™ 


r < 0.430 



more accurately, the 68% and 95% areas under the prob- 
ability curve are covered, respectively, by the following 
intervals of the parameter r: 



€ (0, 0.452) and r G (0, 0.843), 



(13) 



so that the point r = (barely) belongs to the 68% 
interval. 

The 3d ML values, n s = 1.052 and r = 0.285, ob- 
tained with CosmoMC are smaller than those found in 
our previous work Q, n s — 1.111 and R = 0.264 (which 
is equivalent to r — 0.550), in exactly the same setting. 
At the same time, the ML value of A s is larger than 
before. Physically, this means that the new analysis in- 
dicates a larger contribution of density perturbations and 
a smaller contribution of gravitational waves. The out- 
come r = 0.285 in (fTTj) is almost a factor of 2 smaller 
than the previous number r — 0.550. Taken for the face 
value, these new evaluations weaken the indications of 
relic gravitational waves in the WMAP7 data from ap- 
proximately a 2a level to approximately a la level. This 
also worsens the prospects of confident detection of relic 
GW by the Planck satellite (more details in Sec lVI[) . Nev- 
ertheless, it is fair to say that some qualitative indications 
of blue spectral indices, and possibly of a large amount 
of gravitational waves, have survived. 

It is unclear to us why the CosmoMC numerical tech- 
nique has led to a factor of 2 grimmer results for r than 
the previous analysis in [6J. It does not seem likely that 
this happened because of a better treatment of noises. 
The increased noises would probably lead only to a larger 
spread of likelihood functions, but this does not seem to 
be the case. The difference is more like a systematic shift 
of ML points for r by a factor of 2 toward smaller values 
of r. More work is needed in order to understand the 
cause of this discrepancy. 



B. Varied background parameters 

To alleviate the worries about fixed background pa- 
rameters, we run the CosmoMC option allowing to find 
the background parameters together with the perturba- 
tion parameters. The information on the background 
parameters is provided by the widely quoted external 
data (not affecting the perturbation parameters). We 
now set free not only the perturbation parameters r, n s , 
ln(10 1( M s ), but also the background parameters Vl^h 2 , 
JVi 2 , r, 9. Together with the WMAP Collaboration 
[l(| , we include in the code the external data on Hq [2l[ , 



BAO [21|,l22j, and SNIa [23j]. In addition, we impose a 
prior on the "age of the Universe" t e (10, 20)Gyr. The 
CMB data are used only up to £ max — 100. We call this 
set of assumptions the Case IV. 

The analysis has led us to the following marginalized 
Id values of the background parameters: 



Q b h 2 



n n9977+ 01906 

U.UZZ I I -0.01962! 

71 2+0.039 

u - ' lz -0.030- T " 



n c h 2 

0.095: 



0.1246 



+0.0236 
-0.0229' 



-0.067 
0.080' 



They are quite close to the values that we used as fixed 
parameters in the Case III. 

As for the perturbation parameters, the 7d ML point 
was found at 

r = 0.313, n s = 1.075, ln(10 10 A s ) = 3.026, (14) 

and rit = 0.075. We also derived the Id mean results (see 
also Table |nj): 

(15) 



1 064+ 055 



ln(10 10 A s ) = 3.0491°°!. 



The Id likelihood functions for n s and r are plotted by 
black curves, also marked by (4), in Fig. [T] The distribu- 
tion for r again, like in the Case III, shows a maximum 
at r = 0.2. The 68% and 95% intervals are covered, 
respectively, by 



r E (0, 0.332) and r e (0, 0.471). 



(16) 



The reported numbers, as well as the shapes of red 
(Case III) and black (Case IV) curves in Fig. [TJ arc 
pretty close to each other. The more concentrated form of 
the likelihood function for r depicted by the black curve, 
as compared with the red curve, is probably the result of 
removal of artificial covariances between the background 
and perturbation parameters implicit in the fixed back- 
ground parameters assumption of the Case III. The Case 
IV is superior to the Case III in that its assumptions are 
more general. And, yet, the results for n s and r turned 
out to be quite similar. The spread of distributions for r 
shown by red and black curves is still considerable, and 
the likelihoods include (barely) the point r = in the 
68% C.L.. Nevertheless, the red and black curves are 
very different from the olive curve in Fig. [T] ( Case I) . 

The results of the Case I are usually interpreted as a 
nondetection of GW with a firm upper bound on r. In 
contrast to the Case I, the results of data analysis per- 
formed along the lines of Case III, and especially Case 
IV, justify our previous conclusion, based now on Cos- 
moMC calculations, that there does exist a hint of pres- 
ence of relic gravitational waves at the level of r ~ 0.2. 
Of course, this is a hint, but not a reliable detection. 
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TABLE II: Results for n s and r in Case III and Case IV 





Maximum likelihood 


1-d likelihood 




n s 


r 


n s 


r (95%C.L.) 


Case III 


1.052 


0.285 


1 064 +l, ' ubB 


r < 0.843 


Case IV 


1.075 


0.313 


i nfu +u,ubb 

l.UD4_ ()53 


r < 0.471 




FIG. 1: One-dimensional likelihood functions for n s (upper panel) and r (lower panel). In both panels, the olive curves, also 
marked by (1), denote the results for Case I, the magenta curves (2) denote the results for Case II, the red curves (3) for Case 
III, and the black curves (4) for Case IV. In both panels, the blue curves, also marked by (5), denote the results for the Case 
V, i.e. for the case with piecewise power-law spectrum. In the upper panel, the blue solid line (5) shows the likelihood function 
for ni 1 ^, whereas the blue dotted line (5) shows the likelihood function for ni 2 '. In the lower panel, the blue line (5) shows the 
likelihood function for r. 



IV. JUDGMENTS ABOUT n s AND r AS 
FUNCTIONS OF £ max 

The assumption of strictly constant spectral indices is 
a theoretical possibility, but not a necessity, and there- 
fore it should be tested by observations. Here, we show 
how the judgments about n s and r depend on the vol- 
ume of the utilized CMB data. Concretely, we adopt the 
assumption of constant n s and nt, Eq.fl^l, but we include 
in the data analysis only the CMB data up to a certain, 
varied in steps, multipole number £ max : 2 < £ < £ max . In 
other words, the varied ! maj[ of this Section replaces the 
fixed largest £ max mentioned in Sec lII Al 

It would be much too time consuming to run the Cos- 
moMC for each step in £ max . Therefore, for the pur- 
pose of this calculation only, we return to the numerical 



techniques of Ref . [6| . We set free the perturbation pa- 
rameters (r, n s , A a ) and fix the background parameters 
at the same best-fit WMAP numbers {Vl b h 2 = 0.02260, 
fl c h 2 = 0.1123, fl A = 0.72 8, r = 0.087) that were used in 
the Case III of Sec lIII Al The perturbation parameters 
are being determined from the TT and TE CMB data 
sets truncated at a given € max . Since the TT uncertain- 
ties are significantly smaller than the TE uncertainties, 
the major role belongs to the TT data. 

The 3d ML values of n s and r as functions of ^ max are 
shown in Fig. [2] The right panel exhibits the relation 
between n s and r which arises when the common vari- 
able ^ max is excluded. The corresponding marginalized 
Id distributions, showing also the 68% uncertainty bars, 
are plotted in Fig. [3] 

It is clearly seen from Fig. [5] and Fig. [3] that when 
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the £ max is sufficiently large, £ max > 350, the value of n s 
approaches n s = 0.96 and r approaches zero. In other 
words, under the adopted assumptions we recover the re- 
sults of Case I and the conclusions of the WMAP "min- 
imal" model [Tol | . Specifically, at £ max = 350 we get the 
3d ML values n s = 0.963 and r = 0.004. 

However, the gradually decreasing £ max leads to the 
gradually increasing n s and r. As soon as £ max drops to 
f max « 160, the 3d ML points and Id distributions turn 
to distinctly blue spectra n s > 1 and nonzero r. At the 
point £ max = 100 we return to the 3d ML result of @: 
n s = 1.111 and r = 0.550 (equivalent to R = 0.264). The 
sharp drop of r at £ max = 50 is probably the consequence 
of small number of participating data points and strongly 
increased uncertainties, as illustrated in the middle panel 
of Fig. El 

As was shown above, within the same assumptions, 
the CosmoMC technique returns the smaller numbers, 
n s = 1.052 and r = 0.285, at the point f max = 100. This 
was the subject of discussion in Case III of Sec lHI Al 
Nevertheless, the general trend toward blue spectra and 
nonzero r at lower multipoles remains the same, indepen- 
dently of the applied numerical code. 

It is interesting to note that the judgments about n s 
and r that we arrived at satisfy the linear relations: r = 
0.14 + 3.20(n s — 1) for 3d ML parameters (right panel 
in Fig. [2} and r = 0.10 + 3.32(n s - 1) for Id likelihoods 
(right panel in Fig. [3}. The increased scatter of points at 
the right ends of these lines is the reflection of increased 
uncertainties for small values of £ max . 

For every fixed £ max , the uncertainties in n s and r (the 
vertical bars around central points in the left and middle 
panels of Fig. [3} can be interpreted as the "degeneracy" 
between n s and r, in the sense that one and the same 
set of data can be quite successfully described by slightly 
different pairs of n s and r. For example, at € max = 100 
(see the left panel in Fig. 2 of Ref. [6] ) , the "degeneracy" 
is represented by the more or less elliptical 2d-area, with 
the major axis providing a linear relation between n s 
and r. The movement along this line illustrates the fact 
that for a given set of data, and at least for £ max = 100, 
a slightly larger r requires a slightly larger n s , and a 
slightly smaller r requires a slightly smaller n s . In terms 
of the left and middle panels of Fig. [3J this would be 
the movement along the vertical bars around the central 
point at a given £ max . 

It is important to note that the inclined lines in the 
right panels of Fig. [2] and Fig. [3J have a different mean- 
ing. They are the result of movement in the horizontal, 
rather than vertical, direction, that is, they represent the 
central values of n s and r for different sets of data char- 
acterized by different limiting £ max . If the hypothesis of 
a strictly constant n s (and a fixed r) were true in the 
entire interval of accessible multipoles, we would expect 
the inclined lines in the right panels of Fig. [2] and Fig. [3] 
to degenerate to a point, surrounded by some uncertain- 
ties. But this did not happen. We interpret this fact as 
a hint of a genuine, even if a very simple, dependence of 



n s on spatial scale. 

V. PIECEWISE POWER-LAW SPECTRA 

The findings described in previous Sections support 
the proposition that it is unwise to postulate one and the 
same power-law spectrum of primordial perturbations for 
all wavelengths covered by the existing data. It is risky 
to do this independently of the issue of relic gravitational 
waves. But what is more important for us, we have shown 
that if this postulate is imposed, any signs of GW in 
CMB disappear, see Sec|TV] and Case I in Sec lII Al On 
the other hand, it is also true that it is easier to demon- 
strate the drawbacks of a given hypothesis than propose 
a better one. In the absence of a firm guidance, we will 
explore the hypothesis which allows the primordial spec- 
trum to consist of two power-law pieces. This means that 
we adopt the GW and DP power spectra in the form of 
Eqs.([3}-(jB}. We want to show that this hypothesis is very 
much consistent with the indications of GW at lower mul- 
tipoles, as found in Sec lIIH and with the preference for a 
blue spectrum at longer wavelengths and a red spectrum 
at shorter wavelengths, as was discussed in Sec lII Bl and 
Sec [TV] 

To be as close as possible to the already performed cal- 
culations, we use fci = O.OlMpc- 1 in Eqs.©-©. This 
choice of fci corresponds approximately to the multipolc 
£ = 100 (24|. Certainly, there is no reason for the piece- 
wise spectrum to be discontinuous at the wave number 
kt, and Eqs. 13}-© take care of this. In general, the pa- 
rameters and A s ^ are not equal to each other, they 
are linked by the relation shown in Eqs. ([I}-©. [28| 

We again use the CosmoMC sampler. The objective 
is to build the likelihood function in the 8-dimensional 
parameter space, consisting of four background param- 
eters (£lbh 2 , Sl c h 2 , t, 6) and four perturbation parame- 
ters (Lo^lO 10 ^), n ( s \ nf\ r). The parameter Af ] is 

expressible in terms of A s , while the spectral indices 
n s and n t in both parts of the spectrum are related by 
nj = n s — 1. The parameters ni 1 ' and are not nec- 
essarily different; it is the maximum likelihood analysis 
of the data that will tell us their preferred values. Obvi- 
ously, if it is postulated that n s = or if it happens 
that n s = in a particular realization of random tri- 
als, then we simply return to the outcomes of the already 
considered Case I in Sec lII Al We call the introduced set 
of assumptions the Case V. 

We again take into account all the WMAP7 data for 
TT, TE, EE and BB correlation functions in the interval 
2 < I < Ux, where £ max = 1200 for TT, 800 for TE, 
and 23 for EE and BB. The 8d maximum likelihood 
values of the parameters r, rti 1 ', in the Case Fare 
found to be 

r = 0.113, = 1.067, = 0.936. (17) 
The marginalized Id likelihood functions for these three 
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TABLE III: Results for n' 1 ', ni 2) and r in Case V 





Maximum likelihood 


1-d likelihood 






ni 


r 






r (95%C.L.) 


Case V 


1.067 


0.936 


0.113 


l-095ir o b 61 


0.958^ 


r < 0.720 




r = 0.14 + 3.20 { n -1 ) 



50 100 150 200 250 300 350 50 100 150 200 250 300 350 095 1 00 1 05 1 ' 10 115 

L_ I n . 



FIG. 2: The maximum likelihood values of n s and r as functions of t 



parameters are plotted in Fig. [T] by blue curves, also 
marked by symbol (5). The Id results, including the 
68% uncertainties, can be summarized as follows (see also 
Table HH, 



1 nQK+0.062 (2) 

i.uyo_ 061 , n s 



0.958 



0.032 
-0.032; 



r G (0, 0.383118) 



Examining the numerical results $T7\ . (|18|) and blue 
curves (5) in Fig. [T]we can conclude the following. First, 

the values for ris > 1 and < 1 confirm the expecta- 
tion that the preferred shapes of primordial spectra are 
blue at long wavelengths and red at short wavelengths. 
The larger (blue) index ni minus its la does not over- 

f 2) 

lap with the smaller (red) index n s plus its la. The 
value of the (blue) index n^p is about the same at the 
values for n s found in the Case III and in the Case IV 
for the interval 2 < I < 100. And the value of the (red) 

(2) 

index n s — 0.958 is about the same (slightly less) than 
n s = 0.969 found in [6] for the interval 101 < £ < 220. 

The general shape of the likelihood function for r (blue 
curve (5) in the lower panel of Fig. [1} is quite similar to 
the distributions (3) and (4) , plus the expected flattening 
of the likelihood function (5) at small values of r. The 
flattening is expected because there is still a large area in 
the parameter space where ni is equal to and both 



of them are less than 1, see the likelihoods (5) in the 
upper panel of Fig. [1] As we know from the analysis of 
the Case /all these options return very small values of the 
parameter r. This is why the probability of small values 
of r has increased in comparison with the "clean" cases, 
Case III and Case IV, and the likelihood function (5) has 
flattened and became consistent with the hypothesis of 
no gravitational waves, r = 0. 

The hypothesis of a piecewise spectrum (Case V) is 
pretty much in agreement with WMAP7 data and with 
all other findings discussed here. But it also illustrates 
how careful one should be in judgments about the ab- 
sence or presence of gravitational waves in CMB data. 
Hopefully, the forthcoming data of better quality will al- 
low us to make more decisive conclusions. 



VI. PROJECTIONS ON THE PLANCK 
MISSION 

The CosmoMC calculations have tempered our previ- 
ous evaluations @ of relic GW in WMAP7 data. The 
indications became weaker, but we do not think they 
disappeared altogether. The projections on the possible 
Planck findings are anticipated to become worse, but it 
is important to see what they are now. The GW sig- 
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FIG. 3: The marginalized One-dimensional distributions (constraints) of n a and r as functions of £ max - The vertical bars denote 
the uncertainties at 68% C.L. 




FIG. 4: The Planck's signal-to-noise ratio S/N for expected GW signals considered in the text. 



nal may have declined in strength, but the good news 
is that Planck is definitely expected to operate for 28 
months rather than for 14 months. The longer duration 
of observations will partially compensate the suspected 
weakening of the GW signal. 



Following |6j , we use the signal-to-noise ratio 
~N = Ar' 



(19) 



where the uncertainty Ar depends on statistics and vari- 
ous instrumental and environmental noises. All our input 
assumptions about Planck's instrumental noises, number 
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and specification of frequency channels, foreground mod- 
els and residual contamination, sky coverage, etc., are 
exactly the same as in our previous papers [6j. In par- 
ticular, the ability, ranging from excellent to none, of 
removing the foreground contamination is characterized 
by the parameter cr fg = 0.01, 0.1, 1. We also work with 
the pessimistic case, in which cr fg = 1 and the nominal in- 
strumental noise in the BB information channel at each 
frequency is increased by a factor of 4. 

The uncertainty Ar is evaluated with the help of the 
Fisher matrix formalism (26| . Since all three perturba- 
tion parameters r, n s and A s are supposed to be deter- 
mined from the same set of Planck's data at 2 < £ < 100, 
we calculate the rr-element of the inverse Fisher matrix, 
Ar = y/ (F^ 1 )^-. All information channels, i.e. TT, 
TE, EE and BB correlation functions, are taken into 
account. The data from the frequency channels at 100, 
143, and 217 GHz are supposed to be used, and the obser- 
vation time of the Planck mission is taken as 28 months. 

We illustrate in Fig. 2] the results of the revised calcu- 
lation of S/N as a function of r. The two benchmark val- 
ues of r derived from the CosmoMC analysis are shown: 
r = 0.285, which is the ML value of r found in the Case 
III, and r = 0.2, which is the Id peak value of r, which 
is found in the Case III and Case IV. Assuming that 
r = 0.285 or r = 0.2 are fair representations of the real- 
ity, the prospects of discovering relic gravitational waves 
with Planck are still encouraging. 

If the true value of the parameter r is r = 0.285, 
the S/N becomes S/N = 7.43, 6.81, 5.43 for a f s = 
0.01, 0.1, 1, respectively. Even in the pessimistic case, 
the detection looks quite confident, because the expected 
level is S/N = 3.27. On the other hand, if r = 0.2 
is the true value of r, the S/N diminishes to S/N — 
5.87, 5.19, 4.05 for cr fg = 0.01, 0.1, 1, respectively. In 
the pessimistic case, the prospects of detection drop to 
S/N = 2.45. 

To conclude, the CosmoMC-revised values of r make 
the forecasts for the Planck mission worse than previously 
evaluated in (aj. It looks like some luck will be needed 



in the context of the foreground removal. Nevertheless, 
even in the pessimistic scenario, the S/N remains at the 
interesting level S/N > 2. It is also necessary to remem- 
ber that the search specifically for the B-mode of CMB 
polarization by the ground-based and suborbital experi- 
ments (see, for example [13]) provides an important extra 
avenue for the detection of relic gravitational waves. 



VII. CONCLUSIONS 

We have reanalyzed the WMAP7 data with the help of 
CosmoMC package and have shown that the (marginal) 
indications of relic gravitational waves are still present. 
It is vitally important not to overlook relic GW in the 
forthcoming data of better quality. The GW signal is 
weak and its discovery can be done, realistically, only 
by parametric methods. Therefore, a correct theoretical 
model and adequate data analysis techniques are espe- 
cially crucial. We have stressed the importance of look- 
ing for GW in the lower-^ interval of multipoles and the 
dangers of unwarranted assumptions about density per- 
turbations. The imminent release of the results of Planck 
observations will hopefully confirm our expectations. 
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